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ABSTRACT 

Ground-based wide-field surveys for nearby transiting gas giants are yielding far fewer true plan- 
ets than astrophysical false positives, of which some are difficult to reject. Recent experience has 
highlighted the need for careful analysis to eliminate astronomical systems where light from a faint 
eclipsing binary is blended with that from a bright star. During the course of the Trans-atlantic 
Exoplanet Survey, we identified a system presenting a transit-like periodic signal. We obtained the 
proper motion and infrared color of this target (GSC 03885-00829) from publicly available catalogs, 
which suggested this star is an F dwarf, supporting our transit hypothesis. This spectral classification 
was confirmed using spectroscopic observations from which we determined the stellar radial velocity. 
The star did not exhibit any signs of a stellar mass companion. However, subsequent multi-color pho- 
tometry displayed a color-dependent transit depth, indicating that a blend was the likely source of the 
eclipse. We successfully modeled our initial photometric observations of GSC 03885-00829 as the light 
from a K dwarf binary system superimposed on the light from a late F dwarf star. High-dispersion 
spectroscopy confirmed the presence of light from a cool stellar photosphere in the spectrum of this 
system. With this candidate, we demonstrate both the difficulty in identifying certain types of false 
positives in a list of candidate transiting planets and our procedure for rejecting these imposters, 
which may be useful to other groups performing wide-field transit surveys. 

Subject headings: stars: binaries: eclipsing — stars: planetary systems — techniques: photometric — 
techniques: radial velocities — stars: individual: alphanumeric: GSC 03885-00829 



1. INTRODUCTION 

The discovery of the 169 known extrasolar planets 8 
has greatly enhanced our understanding of planetary sys- 
tems. Most of these extrasolar planets have been iden- 
tified from Doppler surveys, which search for the radial 
velocity variation of a star caused by the presence of a gas 
giant. These observations can estimate the period and 
eccentricity of the planetary orbit, and provide a lower 
limit for the planetary mass relative to that of the star. 
Assuming that the stellar mass can be precisely deter- 
mined, perhaps by comparing t he stellar spec tra to stel- 
lar atmosphere models (such as IKuruc z 1993|, the mini- 
mum mass of the planet can be derived. The diversity of 
these planetary systems has drastically altered our the- 
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oretical appreciation of their morphology and evolution, 
including the existence of "hot Jupiters" , Jupiter-sized 
planets with periods of a few days that experience high 
insolation from the nearby star. However, less is known 
about the nature of the planets themselves, except for the 
cases when a hot Jupiter was observed to pass in front 
of a dwarf st ar. By observ ing such a transit (as was first 
suggested bv lStruvelll952j) . we can estimate the radius 
of the transiting planet from the fraction of starlight the 
planet blocks during the transit and the radius of the star 
itself, which must be measured from stellar model fits to 
spectra of the star. For a transit to occur, the plane- 
tary orbital inclination must be approximately 90°; this 
implies that the Doppler limit for the planetary mass 
must be close to the actual mass of the planet. We 
now have estimates for the radii and masses of 9 plan- 
ets from transit observ ations similar to tho s e first suc- 
cessfully performe d by iCharbonneau et aD l)2000|) and 
IHenrv et alJ J2000). These estimates provide constraints 
for the competing theories of planetary formation. For 
example, the recently discovered transiting "hot Saturn" 
HD 149026 b has a mass and radius that imply the planet 
has a large core (approxim ately 70Mcp: ISato et alJl2005t 
ICharbonneau et alJ 12006^ . It is hypothes ized that this 
planet must have fo rmed via core accretion ( Pollac Ul984t 
IPollack et al.lll996l). r ather than through gravitational 
instability (B osall997jh However, not every observation 
of a transiting planet can be explained by current models 
of close-in giant planets experiencing high stellar insola- 
tion. The transiting planet HD 209458 b has a radius 
larger than other planets of its mass, and larger than the 
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predicted ra dii fseelLaueh in et alJl2005l and references 
therein, and lDeming et al! 120051) . 

Several wide-field photometric surveys with the goal 
of identifying nearby transiting planets are currently ac- 
tive. Our Trans-atlantic Exoplanet Survey 9 (TrES) is 
a network of three 10 cm telescopes: Sleuth (located 
at Palomar Observ atory, California), PS ST (Lowell Ob- 
servatory, Arizon a: iDunham et aTll2004D . and STARE 10 
(Tenerife, Spain; lAlonso et al]l2004b|K The TrES cam- 
paign, together wi th other wide-fiel d surveys such as 
the HAT network jBakos et alJ 12002]) and Super WASP 
l)Street et al.ll2003l) . are monitoring thousands of nearby 
bright stars (9 < V < 13). We hope to find recurring 
eclipses with the short period and small amplitude cor- 
responding to a transiting hot Jupiter. The brightness of 
the target stars facilitates both the photometric precision 
and the follow-up of any identified transiting planets us- 
ing space-borne telescopes. Examples of detailed follow- 
up observations include the measurement of several 
chemical abundances in the atmosphere of HP 209458 b 
ijCharbonneau et al.l2002T:lVidal-Madiar et al.l2003j) and 
the first direct detections of emitted planetary radiation 
IjCharbonneau et al]l2005l iDeming et al.ll2005l) . 

Many astrophysical systems exist that mimic the light 
curve of a transiting hot Jupiter. Due to the mass-radius 
degeneracy for bodies with masses between 0.001 and 
0.1 M©, the depth of an eclipse of a solar-type star by a 
hot Jupiter, a brown dwarf or an M dwarf star will be 
about 1% in each case, despite the large range in mass. 
Grazing incidence eclipsing binaries may also exhibit 
comparable eclipse depths. For wide-field ground-based 
surveys, the frequency of the false positives is greater 
than the frequency of detection of true transiting plan- 
ets, by at least an order of magnitude. As part of TrES, 
we typically identify 10-20 of these transit-like photo- 
metric signals out of 15,000-25,000 stars (10 < V < 15) 
in e ach 6° x 6° target field of view (b ~ 15°) (see, 
e.g.. IDunham et aLll2004|) . and similar yields should be 
expected from other wide-field ground-based surveys. 
This number (which is dependent on the density of star 
counts, and hence the Galactic latitude) is consistent 
with theoretical predictions. For example. iBrownl l)2003f) 
predicts that for every 25,000 stars observed, we will find 
10 false positives and only one true transiting planet. 
This assumes that we must observe 3 transit events for 
each candidate, and is dependent on the visibility of tran- 
sits throughout the observation run. The low yield of 
planets necessitates a rigorous routine of follow-up ob- 
servations and detailed analysis to eliminate all possible 
alternatives to the planet hypothesis. One straightfor- 
ward method to reduce the number of false positives is to 
obtain multi-epoch spectroscopy of each candidate and 
measure radial velocities with a precision of ~ 1 kms -1 . 
From this we can identify targets with compan ions of 
stellar, rather than planetary, mass (see, e.g.. iLathaml 
120031 iCharbonneau et al.ll2004|) . We can also estimate 
the luminosity class of the target star to single out and 
reject giants. 

The blending of the light from an edge-on binary sys- 
tem with that from a third star can also be mistaken 
for a transit. IBrownl l)2003l) predicts that half of the 

9 http: //www. astro . caltech . edu/~f tod/tres/ 

10 http : //www . hao . ucar . edu/public/ research/stare/stare . html 



false positives from a typical wide-field survey will be 
of this type; the other half will be grazing eclipsing bi- 
naries. Blends can be much more difficult to identify. 
The faintness of the binary compared to the third star 
can prevent the detection of the radial velocity varia- 
tions of the binary. These variations are also masked 
by the rotationally broadened spectral lines of the ro- 
tationally synchronized binary stars. However, if the 
binary has a significant difference in effective temper- 
ature from that of the third star, the eclipse depths 
should display a strong color dependence, unlike the 
color-independent transits of a (dark) planet across a 
single star. Nevertheless, there has been recent expe- 
rience of blends with color-in dependent eclipse d epths. 
In the case of OG LE-TR-33 (ITorres et al.l l2004bf) and 
GSC 01944-02289 l|Mandushev et alJl2005l) both candi- 
dates (a suspected planet and brown dwarf, respectively) 
showed color-independent eclipse depths, and yet were 
subsequently discovered to be blended systems. Evidence 
for the presence of an eclipsing binary was found from a 
careful analysis of the spectral line shapes, prompting the 
authors to compare the photometric data to simulations 
of blends. OGLE-TR-33 was shown to be a hierarchical 
triple consisting of a bright F6 dwarf and an F4+(K7- 
M0) binary. The blend model for GSC 01944-02289 com- 
prises an F5 primary and a G0+M3 binary. In both 
cases, the similarity in color between the primary star 
and the brightest member of the binary explains the con- 
stant eclipse depth at different wavelengths. The high 
occurrence of such false positives and the difficulty in re- 
jecting them requires a detailed study of candidates be- 
fore any ann ouncement is made, as was done in the case 
of TrES-1 l|Alonso et a,l]l2004afl and OGLE-TR-56b 
(|Torres et al.ll2004aj) . 

Here we discuss a promising candidate, 
GSC 03885-00829, from one of our target fields. 
Initial photometric (J3 an d spectroscopic ([|2Jl mon- 
itoring of this candidate strongly suggested that we 
were observing a Saturn-sized companion transiting a 
solar type star every 1.441 days with a transit depth of 
approximately 6mmag. However, follow-up photometry 
(2J displayed a slight color dependence as might be 
caused by a blend, and we were able to model our 
photometry using simulations of blended eclipsing 
binaries (fJSJl. Our best fit model consists of a bright F 
dwarf and a K dwarf binary, and we were able to identify 
the presence of light from the binary in the spectrum 
of GSC 03885-00829 (©. The faintness of the binary 
system prevents us from identifying the presence of 
asymm etric spectral lines, as was do ne by ITorres et alJ 
(2004bJ and lMandushev et all (J2005). In this case, only 
multi-color observations provided us with the necessary 
evidence to call into question the planetary nature of 
this candidate. 

2. TRES TELESCOPE OBSERVATIONS 

In March 2004, we commenced observations of a 6° x 6° 
target field in Draco. The field is centered on our 
V = 4.8 guide star HD 151613 (a = 16 h 45 m 17:82, 
5 = +56°46'54"7 J2000). Between UT 2004 March 29 
and June 22, we observed this field nightly with Sleuth 
at Palomar Observatory (California), and with PSST 
at Lowell Observatory (Arizona). STARE, in Tenerife 
(Spain), did not observe this field as it was undergoing 
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repairs at this time. A total of 15854 photometric expo- 
sures of 90 s each were obtained through either a Sloan r 
filter (Sleuth) or a Kron-Cousins R filter (PSST). 

We bias-subtracted and flat -fielded the images of our 
target field once the data were transferred from the ob- 
servatory computers. We performed the calibration of 
the Sleuth data using customized 1DL routines; we cali- 
brated the PSST data using the zerocombine, cc dproc, 
and f latcombine tasks in the 1RAF 11 package dTodvl 
49931). We reduced the Sleuth and PSST photomet- 
ric data separately as follows using our difference im- 
age a nalysis (D1A) pipeline ( describ ed inlDunham et all 
2004 and based in part uDon lAlardll2000fl . 

We created our reference image for the field from im- 
ages obtained at low air mass on a photometric night 
during dark time. We obtained a standard list of stars 
from this image using profi le-fitting (PSF) ph otometry 
fDAOPHOT n/ALLSTAR: IStetsorJl98lilT99l . We cal- 
culated the equatorial coordinates (a, 6) of these stars 
by matchi ng a subset with the stars listed in the Tycho- 
2 Catalog (H og et aLll2000(K and then spatially interpo- 
lated all of the science images so that the star coordinates 
from each image matched those from our standard star 
list. 

We produced the master image for the Sleuth data set 
by combining f 9 of our best-quality interpolated images; 
we combined 17 images to create the corresponding PSST 
master image. We subtracted each interpolated image 
from this master image. We used aperture photometry 
on the resultant difference images (using the centroids 
derived for the standard star list) to estimate the flux of 
each star in each image. We produced time series con- 
sisting of the differences between the magnitude of a star 
in the reference image and the magnitude of that star in 
each target image in turn. We decorrelated these light 
curves as follows to remove systematic effects typical of 
wide-field surveys, such as those caused by changing at- 
mospheric conditions throughout the night. We listed 
the stars in order of brightness, and divided the list into 
batches of 500. We computed the least-squares fit to the 
light curve of a given star from a linear combination of 
the other light curves in that batch. We then subtracted 
this least-squares fit from the ligh t curve of the s tar. 

In our previ ous stu dies llAlonso et alJ 12004a: 
iMandushev et all 120051: iCreevev et alJ I2005|L we 
presented separate light curves from one of the TrES 
telescopes. For this field, we combined the two data sets. 
For a given star on the Sleuth standard star list, we 
calculated the angular distances between that star and 
the PSST standard stars. We matched the Sleuth star 
with a PSST star if the angular distance was less than 
5" (0.5 pixels). Due to the difference in the selected 
filter and the telescope pointing between the two sites, 
some stars were unique to a given standard star list, and 
no match was found. We appended the time series for 
each matched PSST star to the corresponding Sleuth 
time series, and reordered the combined time series 
chronologically. The data for the unmatched PSST stars 
were simply added to the resultant data set. 

11 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



In order to reduce the computational intensity of our 
transit search, we averaged the combined time series 
in 9 minute wide bins to obtain 2996 binned observa- 
tions. Since central transits should last 3 hours, this 
did not significantly sacrifice temporal resolution of po- 
tential transit events. For ~10,000 stars, the rms scat- 
ter of the binned data was below 0.04 mag . We per- 
formed a search of the time series of these stars us- 
ing th e box-fitting least squares transit -search algorithm 
(BLS; IKovacs. Zucker. fe Mazehll2002j) to identify peri- 
odic transit events. The BLS algorit hm calculates a Sig- 
nal Detection Efficiency (SDE; see IKovacs et al.l l2002fl 
for each candidate, which denotes how significant the de- 
tection is. We identified candidates based on this SDE, 
followed by a visual inspection. 

Many of the candidates identified for this field show 
V-shaped eclipses or a possible ellipsoidal variability, or 
have large depths, making it likely that they are not tran- 
siting planets, but rather eclipsing binaries. However, we 
promptly identified a promising candidate. When the 
data for this star were folded with the photometric or- 
bital period of 1.441 days (calculated using the BLS al- 
gorithm), the resultant light curve (Figure^ displayed a 
shallow and flat-bottomed transit, and no noticeable el- 
lipsoidal variability out of transit. The depth (~6 mmag) 
and duration (1.4 hours) of the occultation are consistent 
with a Saturn-sized planet transiting a solar type star. 
The SDE for this candidate (~ 20) was high relative to 
that calculated for typical TrES candidates (~ 10-15). 

Supporting data from online catalogs provided further 
evidence of the planetary nature of this eclipse. Us- 
ing the SIMBAD 12 database, we identified our candi- 
date as the star GSC 03885-00829 (s ee Table H The 
infrared colors (2MASS J-K,J-H; iCutri et alJl2003]) 
and optical (B — V) colors of this star are near-solar, 
roughly consistent with the stellar parameters inferred 
from our transit observations. This star displays sig- 
nificant proper motion (26 m as/year from the USNO-B 
Catalog; iMonet et al.l 12003). suggesting it is a nearby 
dwarf. As a first check of the possibility of contamina- 
tion of light from a nearby star, we verified that there 
is no bright star visible on the Digitized Sky Survey 13 
(DSS) images within our aperture radius (< 30"). 

With due enthusiasm, we proceeded to obtain follow- 
up observations of this exciting candidate, with the goal 
of rejecting the possibility that this was not a transiting 
planet. 

3. SPECTROSCOPIC FOLLOW-UP 

We confirmed that GSC 03885-00829 was an iso- 
lated dwarf star by spectroscopically monitoring this 
candidate. We observed GSC 03885-00829, together 
with other candidates from this TrES field, with the 
Harvard-Smithsonian Center for A strophysics (CfA) 
Digital Speedometer l)Lathaml Il992h , operated on the 
1.5 m Tillinghast reflector at the F. L. Whipple Obser- 



12 http : //simbad . harvard . edu /| 

13 The Digitized Sky Survey (http://archive.stsci.edu/dss/ I 
was produced at the Space Telescope Science Institute under U.S. 
Government grant NAG W-2166. The images of these surveys are 
based on photographic data obtained using the Oschin Schmidt 
Telescope on Palomar Mountain and the UK Schmidt Telescope. 
The plates were processed into the present compressed digital form 
with the permission of these institutions. 



4 



O 'Donovan et al. 



vatory (FLWO) on Mt. Hopkins, Arizona. The spec- 
tral coverage was 45 A centered on 5187 A at a resolving 
power of A/AA « 35,000 (a resolution of S^kms" 1 ). 
We observed this particular target on UT 2005 May 18, 
May 20 and May 21, at an orbital phase 0.52, 0.86, and 
0.50, respectively (calculated using the orbital ephemeris 
of the planet; see Tabled). 

Radial velocities were obtained by cross-correlation us- 
ing templates chosen from a library of synthetic spectra 
computed for us by J. Morse and based on the model 
atmospheres of R. L. Kurucz (J. Morse & R. L. Kurucz, 
2004, private communication). The typical precision of 
a single velocity measurement is 0.5 km s . We mea- 
sured the radial velocity to be constant (—38.48 km s _1 
with an rms of 0.28 km s" 1 ) within our errors. These 
measurements indicate that the target star is not grav- 
itationally bound to a massive stellar companion. Var- 
ious stellar parameters were estimated, again by cross- 
correlating these spectra against a grid of templates from 
our spectral library, seeking the best match. Assuming 
a solar metallicity, we estimated the effective tempera- 
ture to be T e ff = 6150 K and the surface gravity to be 
log g w 4.4, which suggested this was a late F dwarf star, 
consistent with the proper motion and photometric col- 
ors, and with our transit hypothesis. The formal stellar 
rotation we derived (usini rj lkms -1 ) is actually below 
our spectral resolution. The surface gravity suggests the 
star is unevolved. This constraint is important, as this 
star lies within the range of effective temperatures for 
which an ambiguity exists as to the corresponding mass 
of the star while on the main sequence, depending on 
the degree of evolution. An illustration of this ambiguity 
is shown in Figure where the two locations denoted 
by the open circle and large filled circle have the same 
effective temperature as GSC 03885-00829 but rather dif- 
ferent luminosities. The particular age of the isochronc 
for this figure was selected to show this difference more 
clearly. The fainter (lower) location corresponds to an 
unevolved main-sequence star (logg = 4.34) of mass 
1.15 M©, whereas the brighter location is for a star near 
the end of the hydrogen-burning phase, and has a surface 
gravity of log g — 3.68 and a mass of 1.61 M©. The radii 
differ by a factor of about 2.4. From this example, we see 
that, without a surface gravity constraint, we cannot be 
sure of the radius of our target star, preventing our dis- 
criminating between a planetary and a stellar transiting 
companion. 

4. PHOTOMETRIC FOLLOW-UP 

Although we had found no evidence of a stellar mass 
companion from our radial velocity measurements, the 
possibility remained that the observed transits were 
in fact the eclipses of a faint binary whose light was 
blended with that from the bright F dwarf due to the 
large pixel sizes of our detectors. Follow-up observa- 
tions of higher angular resolution might resolve such a 
blended system. A possible wavelength dependence of 
the eclipse depth from multi-color observations would 
also provide evidence of a blend. We organized a follow- 
up photometric campaign to observe multiple transits of 
GSC 03885-00829 using D > 20" (0.50m) telescopes and 
through several different filters. 

High precision photometry of GSC 03885-00829 was 
made on UT 2005 June 8 using the 1.2 m FLWO tele- 



scope (Arizona). For this clear, photometric night, we 
used MiniCam, a two CCD mosaic, each array being 
2048 x 4608 pixels. Observations were binned 2x2 for a 
faster duty cycle. A Sloan g filter was used, with an expo- 
sure time of 30 s and a corresponding duty cyle of ~50.5 s. 
The telescope was de-focused to a FWHM of ~15 pixels 
(~9") to allow greater photon counts per given exposure 
without saturation. Spreading the star image also serves 
to reduce pixel-to-pixel variations that may not be com- 
pletely removed by flat-fielding. A total of 561 photo- 
metric measurements of the field were made over a total 
of 7.875 hours. The differential light curve of our target 
star was obtained using aperture photometry of this star 
and one of the observed reference stars. We corrected for 
the effect of differential extinction on the time series by 
fitting the out of eclipse data. The resultant light curve 
was then converted to flux units (see Figure 0}. These 
follow-up photometric observations were made over a 
year after our TrES images. We used this separation in 
time to obtain a more accurate photometric ephemeris 
for our candidate, T C (HJD) = 2453529.833 + 1.44122 x E 
(see Tabled. 

GSC 03885-00829 was observed on the night of UT 
2005 June 13 at the 0.82 m IAC-80 telescope at the 
Observatorio del Teide (Tenerife, Spain), using its 
1024x1024 pixels CCD and a Johnson R filter. To 
achieve better photometric precision, a slight defocus was 
applied, so as to image stars with a FWHM of ~4". Ex- 
posures times of 22 s were used, and the readout time 
using a 2x2 binning was ~10s. The images were bias- 
and flat-field corrected, and aperture photometry was 
applied using the package for opti mal aperture photom- 
etry vaphot iJDeeg fc Dovlel 12001(1 . Nine reference stars 
were used to build an essemble reference star. The dis- 
persion of the data points is larger at the end of the night 
as the star was closer to the horizon. We corrected for 
differential extinction; Figure H| shows the derived differ- 
ential light curve. 

We obtained BV(RI)c observations of 
GSC 03885-00829 on UT 2005 June 5 with the Lowell 
Observatory 42" (1.05 m) Hall reflector in combination 
with a 2K x 2K SITc CCD. A total of 223 exposures of 
the program field were accumulated: four in B, 211 in V, 
three in Rc, and five in Ic- We observed 20 photometric 
standards in th e SA107 SA1 08, SA109, SA112, and 
PG1633 fields l(Landoltl 11992(1 in order to calibrate 
the photometry. We obtained the following values 14 
for the standard magnitudes of GSC 03885-00829 (the 
numbers in the brackets show the number of individual 
measurements used to derive the mean magnitudes): 
B = 11.077 ± 0.001 (4), V = 10.465 ± 0.001 (12), 
R c = 10.094 ± 0.003 (3), and I c = 9.722 ± 0.001 (5). 

These initial attempts to observe GSC 03885-00829 
convinced us that the target star displayed a transit- 
like dip and that nearby stars were not variable. This 
reduced the possibility that this signal was caused by a 
chance superposition of a star with an eclipsing binary in 
the large pixel scale of the TrES detectors. We were also 
able to reproduce our g band observations using a model 
of a Jupiter-sized planet orbiting a near-solar type star, 

14 The errors include the uncertainties in the Landolt photome- 
try and the internal scatter of our photometry, but may not account 
for the total systematic error in the observations. 
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although the ingress and egress appeared to be too long 
in duration. 

It was when we compared the different light curves that 
we realized something was wrong with our assumption 
that these were observations of a transiting planet. The 
depth of the transit appeared to vary with wavelength: 
an eclipse depth of 0.4% in the g band and 0.7% in the R 
band. Prompted by the color dependence of the transit 
depths, we ran various simulations of the eclipses visible 
from GSC 03885-00829 in an attempt to rule out the 
possibility of a blend. 

5. BLEND ANALYSIS 

Light curve fits to the g-band observations were car - 
ried out as described in detail bv iTorres et al.l l|2P04b). 
Briefly, we assumed that the measured brightness of 
GSC 03885-00829 is due to the light of an eclipsing bi- 
nary blended with the light of the F star, so that the 
deep eclipses of the binary are reduced in depth to the 
level that we see. We hypothesized that the three objects 
formed a hierarchical triple system (rather than a by- 
chance alignment), and we to ok their physical prop erties 
from theoretical isochrones bv lGirardi et alJ {2000). The 
mass of the F star (1.15 M©) was constrained from its 
effective temperature derived in ^ with the assumption 
that it is unevolved. This modeling produced a reason- 
ably good match to the measured dip for an edge-on bi- 
nary composed of an early K star eclipsed by a very small 
M star (see Figure 0. In this scenario there is no mea- 
surable secondary eclipse. Although the eclipse depth is 
well reproduced, the predicted duration is slightly shorter 
than that observed. In addition, the brightness expected 
for the K star (~15% of the F star in the optical) is such 
that it would be visible in our spectra. This fit could 
only be improved by increasing the size (and mass) of 
the brightest star to a value inconsistent with the sur- 
face gravity inferred from our spectroscopy. 

The periods for our transit candidates are calculated 
assuming each observed transit is of equal depth. How- 
ever, when a candidate is in fact a blended eclipsing 
binary, the possibility exists that the primary and sec- 
ondary eclipses are similar enough in depth to be con- 
fused by our period-finding technique. The BLS algo- 
rithm may derive a best fit period for such a system that 
is half of the true value. Therefore, we explored blend 
scenarios in which the period is 2 x 1.441 days. Fig- 
ure shows the result of our best fit to the TrES r-band 
data, in which the eclipsing binary is composed of two 
late K stars with masses of 0.67 M© and 0.64 M©, with 
an orbital inclination of 84° to the line of sight. The 
model indicates a slight difference in eclipse depths (see 
Figure HO* ) of about 1 mmag, although this difference is 
only marginally visible in observations themselves. En- 
largements of the two eclipse regions are displayed in the 
lower panels of Figure El According to this fit, the time 
of the center of transit previously derived (see Table is 
in fact a time of secondary eclipse; all of the photomet- 
ric follow-up observations shown in Figure 0] were taken 
during a secondary eclipse. The brighter of the K stars 
(i.e., the primary of the eclipsing binary) has only ~3% 
of the light of the main F star in the optical, and is be- 
low our threshold for spectroscopic detection. Figure 
shows the location of the three stars (the filled circles) in 
the H-R diagram. 



According to the blend model described above, the 
eclipse in the g band (a secondary eclipse) is predicted to 
be shallower than in the r band (0.35% versus 0.7%), as 
we indeed observe, although the measured depth (0.4%) 
is slightly deeper than predicted. We attribute this to 
shortcomings in the isochrones used for the blend mod- 
eling, which are not specifically designed for low-mass 
stars. In particular, missing opacity sources and other 
physical ingredients may affect the t heoretical luminosi- 
ties in the optical (V or q) bands (IBar affe et al.l Il998t 
iDelfosse et alJ 120001 iChabrier et alJ I2005|L whereas the 
red and near infrared magnitudes are presumably more 
reliable. A sign of this is seen perhaps in the pre- 
dicted V — K color for the main F star: the isochrones 
giv e V — K = 1.26 (in the Johnson system as defined 
bv IBessell k, Brett) I1988J) . bluer than the typical color 
of a dwarf of this te mperature, V — K = 1.40 (e.g., 
IBessell k, Brettjll988j) . Other stellar evolution models 
specifically designed for low-mass stars such as those by 
iBaraffe et alJ (1998) appear to give more realistic colors. 
Our F star is predicted to have V— K = 1.41 according to 
those calculations (after transformation of the isochrone 
K ma gnitudes from t he CIT to the Johnson system, fol- 
lowing |L"eggetdlT992T) . very close to the empirical value. 
Unfortunately the IBaraffe et all l)1998|) models are not 
publicly available for the Sloan bands, so we are unable 
to use them in our blend modeling. 

As indicated earlier, infrared magnitudes for 
GSC 03885-00829 are available from the 2MASS 
Catalog. In particular, the measured V — K color of our 
candidate in the Johnson system (Tabled is 1.66 ±0.02 
(using transform ations from the 2MASS system by 
iCarnenterl 120011) . The difference with the color of a 
single F star indicates a significant infrared excess of 
about a quarter of a magnitude 15 . This in itself can be 
taken as evidence of contamination from the light of 
a later-type object, providing further evidence that we 
are dealing with a blend. The computed color of the 
combined light of th e three stars in our model using the 
IBaraffe et all (|1998l) isochrones is V — K = 1.64, which 
agrees quite well with the observations and supports our 
interpretation. Other predicted red and infrared colors 
also match the measured values reasonably well (see 
Figure 0). 

Although it is quite possible that additional fine- 
tuning may improve the small discrepancies noted above 
in the g band and provide a near-perfect fit to all ob- 
servations (to the extent allowed by the accuracy of the 
stellar evolution models and observational uncertainties), 
our goal here has been to show how subtle the signatures 
of a blend can be, and that with careful modeling it is 
possible to demonstrate that they are in fact due to a 
blend scenario, and therefore to reject the candidate. 

6. CONFIRMATION OF BLEND MODEL 

In order to test further our blend hypothesis, we ob- 
served GSC 03885-00829 on UT 15 August 2005 with 
the NIRSPEC infrared spectrograph at the Keck Obser- 
vatory. We observed the target in the if -band spectral 
region centered near 2.293 /im, to search for the presence 

15 Hence, even if we underestimated the error in our V— band 
photometry (see §|1J by an order of magnitude, the resulting prop- 
agation of error would not significantly affect the size of this dis- 
crepancy. 
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of features from the 12CO 2-0 bandhead. Such features 
are very weak for mid-G-type stars, and absent for stars 
with spectral types earlier than GO. Hence, the detection 
of such features would indicate the presence of a cool stel- 
lar photosphere, as predicted by our blend scenario. 

We used a 3-pixel- wide slit, which yields a spectral 
resolution of approximately 25,000. We gathered two 4- 
minutc exposures, between which we nodded along the 
slit by roughly 5". We differenced the two exposures 
to subtract the sky emission and any pixel-dependent 
detector bias. We extracted the order spanning the lo- 
cation of the 12CO 2-0 bandhead by summing over a 
15-pixel-wide band centered on the peak of the instru- 
mental profile. A small number of values in the extracted 
1 -dimensional spectrum were corrupted due to bad pix- 
els in the infrared detector. We replaced these values 
(30 out of 1024 pixels) by interpolation. This region 
contains a large number of telluric methane features. 
We produced a model of these features by modifying 
the electronic version of KPNO/FTS telluric spectrum 
({Livingston fe Wallace! Il991[) for airmass, wavelength- 
solution, and instrumental point-spread function. Di- 
viding our extracted spectrum by this model yields the 
stellar spectrum corrected for telluric absorption. 

In Figure|Sl we plot the resulting spectrum, as well as a 
spectrum of the nearby M3V star GJ 725A (which shows 
very prominent CO features) for comparison. The rela- 
tive intensity between the individual CO features at 2.3 
um does not differ for these two types of stars, although 
the overall amplitude of the features is reduced for hotter 
temperatures. We did not attempt to match the spectral 
type of the NIRSPEC data quantitatively, since our pri- 
mary goal is to exclude the planetary hypothesis. The 
spectrum of GSC 03885-00829 clearly shows the 12CO 2- 
bandhead near 2.293 /um; the relative depth of the band 
head is approximately 6%. The detection of this fea- 
ture confirms the presence of a cool stellar photosphere. 
Also, for the late K dwarfs of our model, the 12CO 2-0 
band head has a depth of approximately 30%, and the 
K dwarfs in our blend model contribute 25% of the K— 
band flux from the system. Hence the expective depth 
of this band head in our spectrum is ~ 7.5%, in rough 
agreement with the observed relative depth. Thus we in- 
terpret this feature as originating in the photospheres of 
the K-stars of the binary. 

The heliocentric Julian Day at mid-exposure was HJD 
2453597.83591, which corresponds to an orbital phase of 
0.09, at which point the expected velocity separation be- 
tween the two K-stars is 89.5 km s -1 . The K-stars are 
likely tidally-locked and hence their spectral features will 
have a usini = 12kms _1 . Since this is smaller than the 
predicted velocity separation at the time of the exposure, 
we might expect to resolve the individual components of 
the observed 12CO 2-0 feature from the two K-stars of 
the binary. And indeed the band head in our spectrum 
shows two clear peaks of similar depth, with a velocity 
separation similar to the predicted separation of the K 
dwarfs. A more careful analysis of this spectrum (using 
for example the t wo-dimensional cross-c orrelation algo- 
rithm TODCOR; iZucker fc MazehllT994|) should recover 
the components more precisely. However, for the pur- 
poses of this paper, it is enough to identify the presence 
of the light from an eclipsing K binary system in our 
spectrum, which rules out the transiting planet hypoth- 



esis. 

7. DISCUSSION 

The difficult task of eliminating any contamination re- 
sulting in a false transit signal is the primary challenge 
currently facing wide field transit surveys. Developing 
this experience will not only enable us to make firm 
detections of transiting Jupiters, but will be extremely 
valuable as we search for Earth-sized planets outside our 
solar system with NASA's Kepler mission. As with the 
ground-based wide-field surveys, the challenge with Ke- 
pler will not only be obtaining the photometric precision 
necessary to observe these minute signals, but rejecting 
all other possible causes of these eclipses as well. 

We have presented here one of our disappointments: a 
candidate that passed all of our initial photometric and 
spectroscopic tests, but was later shown to be a result 
of contaminated light from an eclipsing binary. As such, 
it highlights the difficulty in rejecting all false positives 
from a transit survey. As the components of this eclips- 
ing binary are both faint K dwarfs, the resultant radial 
velocity variations in the light blended with that from 
the nearby F star are not detectable. Also, the resul- 
tant color dependence of the blended light curve is small, 
though observable, and the system displays a color red- 
der than that of an isolated F star. 

The TrES survey readily produces on the order of 10 
candidate transiting planets from each selected field of 
view with 15,000-25,000 stars, the majority of these can- 
didates proving to be false positives. Some of these 
will mimic the expected properties of a planetary system 
quite closely. Based on our experience of frequent blends, 
where the telltale indicators of the stellar components are 
often masked, spectroscopic and even multi-color photo- 
metric follow-up is insufficient to confirm the planetary 
nature of a candidate. An attempt must be made to in- 
terpret the observations as those of a blended eclipsing 
binary, and this interpretation rejected only if the obser- 
vations are not in agreement. Having meticulously ex- 
amined the evidence, we can commit to obtain the radial 
velocity orbit of this firm candidate through high reso- 
lution spectroscopy with a high signal-to-noise ratio. It 
should be emphasized that determining the mass of the 
candidate from such observations is a necessary step in 
identifying a transiting planet. The methods to reject 
astrophysical false positives presented here and by other 
authors cannot be used to confirm the planetary nature 
of a candidate; rather they increase the yield of planets 
from the resource-intensive high-dispersion spectroscopy 
required for such a confirmation. 
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TABLE 1 
Data for GSC 03885-00829 



Parameter Value 



R.A. 


(J2000) 


16 h 52 m 33!7 


Decl. 


(J2000) 


+57°58'27" 


GSC 




03885-00829 


2MASS 




16523368+5758262 


V a 


(mag) 


10.465 ± 0.001 


B - V a 


(mag) 


0.612 ±0.001 


V-R c 


a (mag) 


0.371 ± 0.003 


V-I c 


1 (mag) 


0.743 ± 0.001 


J h 


(mag) 


9.197 ±0.018 


J - H b 


(mag) 


0.337 ±0.023 


J-K a 


b (mag) 


0.440 ± 0.026 


Period c 


(d) 


2.88244 ± 0.00046 




(HJD) 


2453529.833 ± 0.009 


Depth 


(r mag) 


0.006 ± 0.003 



a See J3]for a discussion of errors. 

b From the 2MASS Catalog ICutri et alj|2003T) . 

c The period of the suspected candidate planet 
was l d 44122±0 d 00023. 

d The time of the secondary eclipse of the binary, 
and the central transit time of the candidate. 
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Fig. 1. — Binned TrES r light curve of GSC 03885-00829, folded with the photometric period of 1.441 days computed using the Box-fitting 
Least Squares algorithm. Overlaid is the corresponding box transit model. 
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Fig. 2. — Sample spectrum of GSC 03885-00829 (which includes the Mg I b triplet) and the corresponding cross-correlation function. 




Fig. 3. — IGirardi et al. (2000) model isochrones for solar metallicity and ages ranging from 1 to 4 Gyr. The open circle and large filled 
circle represent two main-sequence stars of the same effective temperature as our candidate GSC 03885-00829 (T c g = 6150 K), but different 
degrees of evolution. They are shown on the 2 Gyr isochrone (heavy line), which maximizes the difference in brightness at this temperature. 
The corresponding radii differ by a factor of 2.4. We must therefore constrain the evolution of our candidate using spectroscopy before we 
can accurately estimate its radius and hence the size of any transiting companion. The three filled circles show the location of each of the 
three members of our final blend model for this candidate on the 2 Gyr isochrone. The bright primary has a mass of 1.15 Mq, and the 
binary component masses are 0.67 Mq and 0.64 Mq. 
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Fig. 4. — Our follow-up photometry of GSC 03885-00829 near the predicted time of transit. The observations were folded using the 
orbital ephemeris of the planet (see Table IT!:) . Each plot is labelled with the corresponding telescope (see text), and the filter bandpass. 
The predicted transit events were observed, but an increase in eclipse depth with increasing wavelength is apparent. 



9.990 



9.995 



10.000 - 



10.005 



10.010 



+ +4+ 




-0.1 



0.0 

Phase 



0.1 



Fig. 5.— Binned g-band observations of GSC 03885-00829 obtained with the FLWO 1.2 m telescope, folded with the BLS period of 1.441 
days. Superimposed is the theoretical light curve from a blend model consisting of the bright F star and a K+M dwarf eclipsing binary 
(see text). Although the depth of the transit is well fit, the predicted duration is slightly shorter than observed. 
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Fig. 6. — TrES r-band photometry of GSC 03885-00829, folded using a period twice that of our candidate planet (i.e., 2 X 1.441 days), 
(a) Best-fit blend model consisting of the bright F star and a pair of eclipsing K dwarfs (sec text). The theoretical curve indicates a small 
difference in depth between the primary and secondary eclipse; (b) Enlargement around the primary eclipse; (c) Enlargement around the 
secondary eclipse. The observed duration of the eclipses is well reproduced by the model. 
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Fig. 7. — Near infrared colors measured for GSC 03885—00829 (dots) compared with the theoretical colors of a single F star (dashed line) 
and the colors from our blend model (solid line; combined light of three stars). The latter is seen to reproduce the observed colors well. 
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Fig. 8.— X-band spectra of GSC 03885-00829 and a nearby M3V star (GJ725A), obtained with the NIRSPEC spectrograph at the 
Keck Observatory. The M dwarf spectrum, with strong CO features, is shown for comparison with our target spectrum, which displays 
the 12CO 2-0 band head at approximately 2.293 (im. The spectrum must therefore include the light of a star with a spectral type later 
than the F dwarf target star, such as the K dwarfs of our blend model. The presence of the CO feature thus rules out the possibility of a 
planetary companion. 



